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ABSTRACT 

We have observed the center of the Local Group dwarf irregular galaxy IC 1613 with WFPC2 
aboard the Hubble Space Telescope in the F439W, F555W, and F814W filters. We analyze 
the resulting color-magnitude diagrams (CMDs) using the main-sequence and giant branch 
luminosity functions and comparisons to theoretical stellar models to derive a preliminary star- 
formation history for this galaxy. We find a dominant old stellar population (aged ~ 7 Gyr), 
identifiable by the strong red giant branch (RGB) and red clump populations. From the (V— I) 
color of the RGB, we estimate a mean metallicity of the intermediate-age stellar population 
[Fe/H] = —1.38 ± 0.31. We confirm a distance of 715 ±40 kpc using the I-magnitude of the 
RGB tip. The main-sequence luminosity function down to I ~ 25 provides evidence for a roughly 
constant SFR of approximately 3.5 x 10~ 4 Mq yr _1 across the WFPC2 field of view (0.22 kpc 2 ) 
during the past 250-350 Myr. Structure in the blue loop luminosity function implies that the 
SFR was ~50% higher 400-900 Myr ago than today. The mean heavy element abundance of 
these young stars is ~ 1/10 solar. The best explanation for a red spur on the main-sequence 
at I ~ 24.7 is the blue horizontal branch component of a very old stellar population at the 
center of IC 1613. We have also imaged a broader area of IC 1613 using the 3.5-meter WIYN 
telescope under excellent seeing conditions. The WIYN color-magnitude diagram reveals a 
prominent sequence of asymptotic giant branch stars and red supergiants that is less prominent 
in the WFPC2 CMD due to the smaller field of view. The AGB-star luminosity function is 
consistent with a period of continuous star formation over at least the age range 2-10 Gyr. We 
present an approximate age-metallicity relation for IC 1613, which appears similar to that of 
the Small Magellanic Cloud. We compare the Hess diagram of IC 1613 to similar data for three 
other Local Group dwarf galaxies, and find that IC 1613 most closely resembles the nearby, 
transition-type dwarf galaxy Pegasus (DDO 216). 
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Subject headings: galaxies: individual (IC 1613) 
stellar content — Local Group 

1. INTRODUCTION 

The study of stellar populations is a way for as- 
tronomers to explore galaxy evolution, providing di- 
rect information on the variations of star-formation 
rate and mean abundance with time, from the 
present-day back to the initial formation of a galaxy. 
The ensemble of Local Group dwarf galaxies, because 
of their diversity of sizes, environments, and chemi- 
cal abundances makes a good laboratory in which to 
investigate the processes which control galaxy evolu- 
tion (e.g., Hodge 1989; Da Costa 1998; Grebel 1998; 
Tolstoy 1999). 

Dwarf galaxies are the most common class of galaxy 
in the Universe, but there are a number of unresolved 
questions of dwarf galaxy evolution; the relation be- 
tween the early- (elliptical/spheroidal) and late- (ir- 
regular) type galaxies, the effects of environment on 
galaxy evolution, and the patterns of chemical enrich- 
ment are unanswered questions (Mateo 1998). Addi- 
tionally, it may be that the ubiquitous faint, blue 
galaxies seen in redshift surveys (e.g. Ellis 1997) are 
the precursors of today's dwarf galaxies. The most 
straightforward way to trace this possible cosmolog- 
ical connection is by a determination of the star- 
formation histories of the nearby dwarf galaxies (Tol- 
stoy 1998). 

Color-magnitude diagrams which include one or 
more main-sequence turnoffs are a powerful tool for 
the determination of detailed star- formation histories 
for galaxies. However, stellar crowding at the level 
of the horizontal branch and below has made this 
task all but impossible for ground-based telescopes. 
Consequently, most studies to date have concentrated 
on the sparse, late-type dwarf galaxies, or the rela- 
tively uncrowded outer fields around dwarf irregulars. 
WFPC2 has made it possible to push stellar pho- 
tometry down below the magnitude of the horizon- 
tal branch throughout the Local Group, opening new 
vistas for the study of stellar populations in dwarf 
irregular galaxies. 

In a Cycle 6 Guaranteed Time Observer program, 
we obtained deep BVI images of IC 1613 (DDO 8), 
one of the charter members of the Local Group (Hub- 
ble 1936), in order to constrain the ages and metallic- 
ities of the stellar populations in this object. In this 
paper we present the WFPC2 VI color-magnitude di- 
agram of IC 1613, complemented by ground-based 
images taken at the WIYN telescope. 1 Section [2] gives 
a summary of IC 1613's properties, and in section ^ 
we describe our observations and measurement pro- 

1 The WIYN Observatory is a joint facility of the University of 
National Optical Astronomy Observatories. 
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cedures. Section |] presents our color-magnitude dia- 
grams of IC 1613, and the stellar populations analysis 
is given in section |5[ We discuss IC 1613's range of 
metal abundance and relation to a sample of Local 
Group dwarfs for which WFPC2 stellar population 
studies have been made in section ||, and summarize 
our results in section [j]. 

2. IC 1613: BACKGROUND 

IC 1613 was discovered by Wolf in 1906, who noted 
an extended grouping of faint nebulosities near the 
star 26 Ceti (Wolf 1907). IC 1613 was first resolved 
into stars, at magnitude 17-18, by Baade (1928), who 
identified it cLS Si Magellanic Cloud-type galaxy and 
compared it to NGC 6822. IC 1613 was Baade's 
(1963) prime example of the eponymous sheet of Pop- 
ulation II red giant stars that seem to exist in all 
Local Group dwarfs (Baade 1963; Sandage 1971b). 

Due to low values of both foreground and internal 
reddening along IC 1613 sightlines, it has played an 
important role in the calibration of the Cepheid vari- 
able Period-Luminosity relation for the determina- 
tion of extragalactic distances (Baade 1935, Sandage 
1971a). IC 1613 is the only Local Group dwarf ir- 
regular galaxy, apart from the Magellanic Clouds, in 
which RR Lyrae type variables have been detected 
(Saha et al. 1992; Mateo 1998). 

IC 1613 occupies a relatively isolated position in 
the Local Group, «400 kpc from M33, «500 kpc 
from M31, and ~700 kpc from the Milky Way. Its 
known neighbors within 400 kpc are the compara- 
bly luminous WLM irregular (DDO 221), at «390 
kpc, the smaller Pegasus dwarf (DDO 216), some 360 
kpc distant, and the tiny (My = —9.7) Pisces dwarf 
(LGS 3), at ~280 kpc distance. This minor quartet of 
galaxies share similar ranges of chemical abundance, 
and a roughly constant long-term star- formation rate. 
They show varying amounts of recent star formation, 
but each have stellar populations that in the mean are 
older than ^ 3 Gyr (Pegasus — Gallagher et al. 1998; 
Pisces— Aparicio, Gallart & Bertelli 1997; WLM— 
Minniti & Zijlstra 1996; IC 1613— Freedman 1988a; 
this paper). Pegasus and Pisces have been classified 
as "transition" galaxies, intermediate between the ir- 
regular and spheroidal classes (Mateo 1998), although 
the significance of this distinction has been called into 
question (Aparicio et al. 1997). The newly discovered 
Andromeda VI (Pegasus dwarf spheroidal) also lies ss 
400 ±40 kpc from IC 1613 (Grebel & Guhathakurta 
1999). 

The basic parameters of IC 1613 are summarized in 
Wisconsin-Madison, Indiana University, Yale University, and the 
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Table 1 
Basic Properties of IC 1613 



Property Value Source 



£, b (2000.0) 


129?8, -60?6 


Mateo 1998 


(m-M)o 


24.27 ±0.10 


Lee et al. 1993 


Eb-v 


0.03 ±0.02 


Mateo 1998 


My 


-15.2 ±0.2 


Hodge 1978 


B-V 


0.71 ±0.1 


Hodge 1978 


M(HI) (10 6 M Q ) 


54 ± 11 


Lake & Skillman 1989 


M(HI)/L B (M /L Q ) 


0.81 


Mateo 1998 


Rcore (arcsec) a 


200: 


Hodge et al. 1991 


L (Ha) (erg s _1 ) 


3.18 x 10 38 


Hunter et al. 1993 


SFR (10~ 3 M /yr) b 


2.9 ± 0.7 


this paper 



a core radius of best fit King model to unresolved light profile. 
b derived from extinction-corrected Ha flux, see text. 



Table 1. The distance is unusually secure due to the 
large number of Cepheid variables and low reddening 
towards IC 1613 (Sandage 1971a, Freedman 1988b); 
the RR Lyrae distance (Saha et al. 1992) and RGB 
tip distance (Freedman 1988a; Lee et al. 1993) are in 
agreement to within the errors. The reddening maps 
of Burstein & Heiles (1982) show E B - V < 0.03 mag 
towards IC 1613; Freedman (1988a) derived Eb-v = 
0.04. Recent work by Schlegel, Finkbeiner & Davis 
(1998) finds a foreground component Eb-v = 0.01. 

We include the scale length of IC 1613's surface 
brightness profile in Table 1, to emphasize the fact 
that the WFPC2 field, ~150" across, images only a 
small portion of the central galaxy. As Hodge et al. 
(1991) emphasize, the tabulated scale length is highly 
uncertain due to the low light levels, small number 
statistics, and the hard-to-define tradeoffs between 
model parameters. 

The Ha luminosity in Table 1 is taken from Hunter, 
Hawley, & Gallagher (1993), adjusted downwards to 
compensate for their assumed distance of 900 kpc to 
IC 1613. The Ha luminosity of a galaxy provides 
a method to estimate the recent (0-10 Myr) star- 
formation rate. Hunter & Gallagher (1986) derive 

M = 7.07 x 10~ 42 L(Ha) M yr -1 , (1) 

with a possible additional factor of ~1.6 to be added 
to account for internal extinction within their sample 
of calibrating galaxies. We find a total star-formation 
rate over the past 10 Myr of 2.9 ±0.7 x 10~ 3 M yr -1 
for IC 1613. 



3. THE DATA 

3.1. Observations & Reductions 

IC 1613 was observed with WFPC2 over 9 orbits 
on 26-27 August, 1997. The PC chip was centered 
97" southwest of the center of the galaxy, at right 
ascension l h 04 m 48!7, declination ±02° 07' 06"2 
(J2000.0). Exposures totalling 10,700 sec each were 
taken through the F555W (WFPC2 V) and F814W 
(WFPC2 I) filters. F439W (WFPC2 B) images to- 
talling 2,600 sec were also obtained. The short, blue 
exposures were intended to provide photometry of 
the younger, more luminous main-sequence and blue 
loop stars. The images are available online at the 
HST Data Archive under the listing for Program ID# 
6865. 

All images were reduced using the standard 
pipeline at the Space Telescope Science Institute 
(Holtzman et al. 1995a), and combined using 
a cosmic-ray cleaning algorithm. The combined 
F555W image is shown in Figure [I]. The field of 
view was oriented to place the center of IC 1613 in 
the WF3 chip; most of the central H I hole (Skill- 
man 1987; Lake & Skillman 1989) is contained within 
the five square-arcminute field. Additionally, the field 
was chosen to avoid H II regions (Hodge, Lee & Gur- 
well 1990) in order to focus attention on the older 
main-sequence stars that cannot be photometered 
from the ground due to crowding. The low stellar 
and gas density of IC 1613 are clearly evident by the 
clarity with which background galaxies can be seen. 
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Figure 1 

can be obtained from 
the following URL: 



http://www.astro.wisc.edu/ ~cole/icl613.figl.eps 



FlG. 1. — WFPC2 F555W image of the center of IC 1613. North and East are as marked; the field of view is roughly 3.'5 across 
the widest corner-to-corner span. The low stellar density of IC 1613 is highlighted by the large number of background galaxies in 
the image. 



3.2. Photometry 

We performed profile-fitting photometry on our re- 
duced, combined images using the DoPHOT routines 
and techniques described in Schechter, Mateo and 
Saha (1993); the code has been modified to handle 
the complexities of the WFPC2 point spread func- 
tion (Saha et al. 1994). Using the DoPHOT star- 
matching routines, we identified 47,263 stars in both 
the F555W and F814W frames; the shallower F439W 
exposures permitted the cross-identification of just 
7,922 stars in the F439W and F555W images. Our 
data have limiting magnitudes for photometry of V 
~ 26.5, I - 25.5, and B ~ 25. 

We transformed our photometry to the Johnson- 
Cousins BVI system by making aperture corrections 
to the 075 aperture recommended by pdoltzman e\ 
al. (1995~b)l , adopting the STMAG system zeropoints 
therein, and converting the resulting values to the 
BVI system using the iterative procedure described 
in Holtzman et al. (1995b")| , with coefficients from 
their Table 10. Figure |2| contains our final WFPC2 
(I, V— I) color-magnitude diagram for IC 1613. Typ- 
ical 1-0" errorbars at 24th magnitude are ay = 0.03, 
o"I = 0.04. The distributions of magnitude vs. error 
are shown in Figure ||. 

Along with magnitudes and standard errors, 
DoPHOT returns for each object a "quality" parame- 
ter, which allows the astronomer to discriminate true 
point sources from unresolved blends, background 
galaxies, and stars contaminated by cosmetic CCD 
blemishes. Because of the steep main-sequence lumi- 
nosity function, the blending criterion most strictly 
winnows the photometric sample at faint magnitudes. 



For this paper, we have retained in our analysis 
only those stars with the highest quality photome- 
try ("Type 1"). The remainder of this paper con- 
siders only the 9,810 stars with superior F555W and 
F814W photometry, and the 1,639 stars with superior 
F439W photometry. Because of the small numbers of 
luminous, blue stars in our chosen field, the F439W 
photometry has been used solely to check our adopted 
reddening values for internal consistency. 

The smoothly rising error curves in Figure || in- 
dicate that color-magnitude diagrams will not be 
systematically distorted by the effects of crowding 
(Dohm-Palmer et al. 1997a, Han et al. 1997). From 
experience with WFPC2 photometry, we estimate 
our 90% and 50% completeness levels in (V,I) fall at 
(23.8, 23.8), and (25.6, 25.2), respectively. Complete- 
ness issues are most problematic for the lower main- 
sequence and for the giant branch below the level of 
the red clump. Completeness is high over most of 
the magnitude range in Figures || and [3L decreasing 
smoothly towards fainter magnitudes; because most 
of our analysis deals with stars brighter than 24th 
magnitude, incompleteness is not expected to dra- 
matically affect our results. 

3.3. Observations with WIYN 

To broaden our areal coverage of IC 1613, we ob- 
tained V and I images centered on the WFPC2 posi- 
tion with the 3.5-meter WIYN telescope at Kitt Peak, 
on 1-3 October, 1997. The best data, comprising 
2x600 seconds in I and 900 seconds in V, were taken 
under photometric conditions with 076 seeing. Ad- 
ditional exposures with ~1" seeing totalling 1500 sec 



V-I (mag) 

Fig. 2. — WFPC2 I, V— I color-magnitude diagram, including only stars for which DoPHOT achieved the best possible pho- 
tometry. 9800 stars are shown. The most striking features of this CMD are the prominent red giant branch and the upper 
main-sequence. 



Figure 3 

can be obtained from 
the following URL: 



http://www.astro.wisc.edu/~cole/ic 1613.fig3.eps 



Fig. 3. — The distributions of magnitude vs. photometric error for our WFPC2 BVI photometry. The mean errors at 24th 
magnitude are av = 0.03, ai = 0.04, as = 0.03. 
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younger counterparts to the RC stars. Crowd- 
ing in Figure [5] renders them indistinguishable 
from the RGB. 

5. The asymptotic giant branch (AGB), extend- 
ing redward from the tip of the RGB, is weakly 
present in Figure [2L but more strongly appar- 
ent in the WIYN CMD. AGB stars indicate an 
intermediate-age population (2-6 Gyr). A hor- 
izontal AGB indicates metallicities of ~ 1/10 
solar; additionally, a more metal-poor compo- 
nent lies along an upward extension of the RGB 
in Figure ||. 



6. 



7. 



in each filter were obtained. The broad (6'8 square — 
~9 times the area of the WFPC2 field) field of view 
and fine image quality of the WIYN imager allow us 
to extend our CMD analysis to luminous, short-lived 
phases of stellar evolution which are typically under- 
represented in the WFPC2 field. A WIYN image of 
IC 1613, with the WFPC2 field overlain for compar- 
ison, is shown in Figure |I[ 

The WIYN images were overscan corrected, 
bias-subtracted, and flat-fielded within IRAF's 
ccdproc routines 2 . We performed photometry us- 
ing DoPHOT, as per the WFPC2 images, obtaining 
"Type 1" photometry for 14,678 VI pairs. We also 
observed standard fields from Landolt (1992) to re- 
duce the data to a standard system for comparison 
to the WFPC2 data. The WIYN color-magnitude 
diagram is shown in Figure [|. Comparison to Fig. 
|H shows good agreement regarding the magnitude of 
the tip of the red giant branch and the colors of the 
red giant branch and main-sequence. Comparison of 
individual stars on the red giant branch finds Iwiyn 

= I\VFPC2 —0.04 ±0.03, (V— I)wiYN = (V-I)wFPC2 

+ 0.04 ±0.04. 

4. COLOR-MAGNITUDE DIAGRAMS 

Figures |2| and [5] show two main branches, and sev- 
eral smaller features, each of which contains stars 
with a range of ages and metallicites: 

1. The red giant branch (RGB), between V— I = 
0.7-1.7 and extending upwards to I = 20.3; the 
strongest feature of the CMD, contains stars 
from 500 Myr to ;>10 Gyr among its vari- 
ous substructures. The location of the RGB 
tip agrees well with the adopted distance to 
IC 1613. 

2. The blue plume, blueward of V— I = 0.2, 
consisting primarily of main-sequence stars 
younger than ~1 Gyr. The majority of blue 
plume stars in Figure |2| are expected to be of 
spectral class B; Figure [| includes some prob- 
able O stars, indicators of very recent (^ 10 
Myr) star formation. 

3. The red clump (RC), centered at I = 23.76, V— I 
= 0.9. This feature contains core-helium burn- 
ing stars with ages between ~ 1-10 Gyr, and 
accounts for a plurality of the stars in Figure ^. 

4. The blue loop (BL) stars, extending vertically 
from the RC to I ~ 22.5 before turning di- 
agonally blueward to become indistinguishable 
from the main sequence at I ~ 21. These are the 

2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



The red supergiants (RSG), extending nearly 
vertically to I = 16.5 at V— I ~ 1.3 in Figure ||. 
Fewer than 16 of these are expected to be fore- 
ground dwarfs (Bahcall & Soneira 1980); their 
possible presence does not compromise the real- 
ity of the IC 1613 RSG sequence. These young 
stars are on their way to an 'onion-skin' nuclear 
burning structure and are younger than a few 
10 7 years. Such young stars should be repre- 
sentative of the most metal-rich population in 
IC 1613. 

The horizontal branch, visible at the base of 
the RC blueward extension at I = 24.3, 
and possibly red extension from the main- 
sequence at I = 24.7. These are the most metal- 
poor, oldest (10-15 Gyr) core-helium burning 
stars in IC 1613, tracers of a population whose 
main-sequence turnoff would be expected to lie 
at I ss 28. The gap in the horizontal branch 
between 0.3 < V— I < 0.6 is consistent with the 
location of the RR Lyrae instability strip in the 
CMD. 



The WFPC2 and WIYN color-magnitude diagrams 
contain stars as young as the dynamical timescale 
for protostellar collapse and potentially as old as the 
Galactic globular clusters. The salient features of 
the CMDs are summarized in Figure ||. The divid- 
ing point between the blue and red helium burners 
(HeB) has been defined here to be the Cepheid insta- 
bility strip, by analogy to the blue and red horizontal 
branches. Each feature is marked according to the 
approximate logarithmic age expected for the given 
evolutionary phase; the prominence of the marker 
symbols provides a rough guide to the strength of the 
components of the CMD. Taken together, the CMDs 
provide an overall picture of the star-formation his- 
tory of IC 1613: it appears to have formed stars over 
wide ranges of age, perhaps continuously; most of the 
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Figure 4 

can be obtained from 
the following URL: 

http://www.astro.wisc.edu/~cole/ic 1613.fig4.eps 



FlG. 4.— 900-second, V band, WIYN image of the main body of IC 1613. North is up, East is to the left, and the WFPC2 field 
of view is overlain. The seeing is 0'.'6, causing the background galaxies to stand out clearly through the stars of IC 1613. 
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V-I (mag) 

FlG. 5. — WIYN I, V— I color-magnitude diagram, including only stars for which DoPHOT achieved the best possible photom- 
etry. 14678 stars are shown. Image crowding has blended the red giant branch with blue loop stars and red supergiants, causing 
an artificially large color spread. The luminous red supergiants and extremely red tail of the asymptotic giant branch are more 
fully sampled by the 6'8 x 6!8 field of view. 
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Fig. 6. — Reference card for the stellar populations at the center of IC 1613 as derived from the morphology of WFPC2 and 
WIYN color-magnitude diagrams. CMD features are associated with populations of various ages. The relative strength of each 
population is given by the shaded circles in the appropriate logarithmic age bins. A general picture of approximately continuous 
star formation emerges, although our ability to resolve "bursty" behavior becomes progressively poorer with age. 
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star-formation may have taken place a few to ten bil- 
lion years ago. In subsequent sections of this paper, 
we examine the implications of CMD morphology for 
the star-formation history of IC 1613 in more detail. 

We take the distance to IC 1613, from Lee et al. 
(1993), to be 715 kpc, corresponding to a distance 
modulus (m-M)o = 24.27. We correct for a small 
amount of interstellar reddening, Eb-v = 0.03 (Ma- 
teo 1998). We derive a value E B - V = 0.02 ±0.02 
from a comparison of the B— V and V— I colors of the 
upper main-sequence between 24 > V > 23, using 
the intrinsic color of the zero age main-sequence as 
a template. We adopt the Galactic reddening law of 
Cardelli, Clayton & Mathis (1989), using R v = 3.1, 
to derive extinction parameters Ay = 0.093, At = 
0.045, and E v _i = 0.048. 

5. THE STELLAR POPULATIONS OF IC 1613 

The mean color of the red giant branch, in a simple 
stellar population, can be used to estimate the metal- 
licity of its stars. Da Costa & Armandroff (1990) 
give a calibration for globular clusters' metallicities 
using the mean, dereddened V— I color of their gi- 
ant branches at an absolute magnitude Mi = —3. At 
the distance and reddening of IC 1613, this trans- 
lates to I = 21.25. From the WFPC2 CMD, we 
use 47 stars to determine V — Iq = 1.29 ± 0.07; 
the much larger sample of 255 stars in the WIYN 
field yields V-I = 1.28 ± 0.10. The Da Costa 
& Armandroff (1990) calibration gives [Fe/H]wFPC2 
= -1.38 ± 0.31, and [Fe/H]\viYN = -1-43 ± 0.45. 
These abundance spreads are not solely due to pho- 
tometric errors, which limit the resolution in [Fe/H] 
to 0.17 dex for WIYN, and 0.13 dex for WFPC2; 
an intrinsic spread in color is present due to ranges 
in age and/or metallicity. Our estimate of [Fe/H] ~ 
— 1.4 is in excellent agreement with the previous as- 
sessment [Fe/H]RGB = —1.3 ±0.8 (Freedman 1988a). 
Although IC 1613 is certainly not a simple, single 
age stellar population, this calculation still provides 
a useful first estimate of the stellar metallicities (e.g., 
Da Costa 1998). 

A real abundance variation is almost certainly 
present in IC 1613. However, the presence of 
intermediate-age (5-10 Gyr) stars, strongly suggested 
by the prominent red clump at I ~ 23.8, compli- 
cates the estimation of metallicity from RGB color. 
In stellar populations that are significantly younger 
than the globular clusters, the giant branches lie to 
the blue of the Da Costa & Armandroff (1990) cali- 
bration sequence, causing the V — Iq method to un- 
derestimate the metallicity by up to as 0.3 dex (see, 
e.g., Da Costa 1998). This is the infamous age- 

3 The luminosity function of IC 1613 is thus nearly identical to 
al. (1997a) found N(m) oc xa° v AA . 



metallicity degeneracy; the inferred metallicity range 
among IC 1613's red giants may be, in part, due to an 
age range. If age and metallicity are correlated (e.g., 
with chemical enrichment occuring over time) , than it 
is possible to underestimate the range of metallicities 
among giant branch stars. Independent abundance 
information, together with detailed modelling of the 
star-formation history, is needed in order to lift the 
age-metallicity degeneracy. 

To test our photometry, we rederived the distance 
to IC 1613 using the tip of the RGB in our WFPC2 
CMD. We followed the metallicity calibration given 
by Lee et al. (1993), and applied a Sobel filter (Sakai, 
Madore & Freedman 1996) to the I-band luminosity 
function. This fixed the tip of the RGB at I = 20.28 
±0.09, which yields a distance modulus (m— M)o = 
24.29 ±0.12 for E B _v = 0.03 ±0.02 and [Fe/H] = 
— 1.38 ±0.31. This is in excellent agreement with the 
adopted value of (m-M) Q = 24.27 ±0.10. 

5.1. Recent Star- Formation History 

The young stellar populations of IC 1613 have 
been studied by Sandage & Katem (1976), Hodge 
(1978), Freedman (1988a,b), Hodge et al. (1991), and 
Georgiev et al. (1999). The picture that emerges 
from these studies is of a galaxy in which the cen- 
tral regions have undergone continuous star forma- 
tion over the past 1-300 Myr, with the youngest stars 
clustered into small associations and possible clus- 
ters (Hodge 1978). While the regions of recent star 
formation are more centrally concentrated than the 
underlying old population, the very youngest stars 
seem to avoid the dead center of IC 1613 (Hodge et 
al. 1991). This partially accounts for the lack of blue 
stars brighter than Mi « -3.5 in the WFPC2 field of 
view. 

Freedman (1988a) found that the main-sequence lu- 
minosity function brighter than My = 21 could be fit 
with a simple power law, similar to other Local Group 
dwarf irregulars, with N(m) oc my 58±0 ' 07 . Freed- 
man attributes a flattening of the luminosity function 
slope at V ~ 21 to crowding in her data. The WFPC2 
V-band luminosity function is plotted in Figure |7[ We 
find N(m) oc my 48±0 ' 09 , marginally consistent with 
the steeper slope from Freedman (1988a) 3 . Statistical 
noise begins to dominate the WFPC2 main-sequence 
luminosity function brighter than V ~ 21.5; incom- 
pleteness due to crowding becomes problematic below 
V rj 24.5. Within these limits, which roughly bracket 
the range of the B stars (3-20 Mq; ages between 
10-350 Myr), the constant slope of the logarithmic 
luminosity function suggests an approximately con- 
stant star-formation rate. We integrate a Salpeter 

that of the Sextans A dwarf irregular, for which Dohm-Palmer et 
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Fig. 7. — Main-sequence V-band WFPC2 luminosity function for all photometry. Incompleteness becomes a serious problem at 
V ~ 24.5. The least-squares fit slope is —0.48 ±0.09; this includes stars as old as ~ 350 Myr. Normalized for a Salpeter IMF this 
implies a mean SFR over that time of 3.5 x 10 -4 Mq yr _1 across the 0.22 kpc 2 WFPC field of view. 



IMF from 0.1-120 Mq to normalize the upper main- 
sequence luminosity function, and find a total star- 
formation rate of 3.5 x 10~ 4 Mq yr -1 across the 
WFPC 2 field. At the distance of IC 1613, this corre- 
sponds to 1.6 x 10 -3 Mq yr _1 kpc -2 . 

For comparison, this is slightly smaller than the 
current star-formation rate in the brightest H II re- 
gion (S10, a.k.a. complex 67), although S10 fills just 
one-tenth the area of the WFPC2 field and therefore 
sustains a much higher intensity of star-formation. 
This behavior is common to dwarf irregular galaxies; 
for example, Dohm-Palmer et al. (1997b) found that 
star formation in Sextans A has occurred in local- 
ized regions ~ 10 2 -10 3 parsecs in size, and that these 
regions last for 100-200 Myr at a time. 

Davidson & Kinman (1982) and Kingsburgh & Bar- 
low (1995) have studied the Wolf-Rayet star DR 1 and 
its surrounding H II region; this object, presumably 
a few million years old, has an abundance between Z 
« 0.004 ([Fe/H] = -0.7, Kingsburgh k Barlow 1995) 
and [O/H] = -1 (Davidson & Kinman 1982). These 
are in reasonable agreement with the H II region neb- 
ular abundances reported by Skillman, Kennicutt, & 
Hodge (1989) of [O/H] « -1.0. Recent B star spectra 
(Kudritzki 1999) provide further support for the case 
that IC 1613 is currently forming stars with nearly 
the same metallicity as the Small Magellanic Cloud 
(Hodge 1989). The region of our present study was 
chosen to avoid known H II regions, and it is of in- 
terest to determine the highest metallicity attained 
during recent (but not current) star formation. 



We can estimate the metallicity of the youngest 
stellar population in the WFPC2 field of view from 
the color extent of the blue loop stars — massive, core- 
helium burning stars — and the slope of the red su- 
pergiant sequence in the color-magnitude diagram. 
The WIYN CMD shows a well-populated RSG se- 
quence, rising to I ~ 16.5; the spread in color of this 
sequence can be due to age as well as metallicity. 
From a comparison to theoretical isochrones (Girardi 
et al. 1999), the majority of the RSGs cannot be as 
metal-rich as Z = 0.004 ([Fe/H] « -0.7); the Z = 
0.001 ([Fe/H] rs —1.3) isochrones provide a better fit 
to the data. 

While the RSG sequence above I « 20 is underpop- 
ulated in the WFPC 2 CMD, the blue loop stars are 
prominent. We used a Monte Carlo code to compare 
the isochrones of Girardi et al. (1999) to this por- 
tion of the IC 1613 Hess diagram. 4 Figure || shows 
the results of this simulation. We have sampled from 
the isochrone set assuming a constant star forma- 
tion rate and Salpeter (1955) initial mass function 
for two metallicities, Z = 0.001 (left panel), and Z 
= 0.004 (right panel). In each panel of Figure || the 
star-formation rate was constant from 0.25 - 1 Gyr. 
The Hess diagram of the WFPC2 data is overplotted 
against the simulations, with isopleths from 2 - 256 
decimag -2 ; each successive isopleth shows twice the 
CMD density of the previous one. 

Figure |8| affords two main points of comparison be- 
tween models and data. Firstly, the red HeB (blue 
loop) stars between I « 23.5 and 22, with V— I 



4 a contour plot of the number of stars per unit magnitude per unit color in the CMD — see Hess (1924); Trumpler & Weaver 
(1953) 
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Figure 8 

can be obtained from 
the following URL: 



http:/ /www. as tro.wisc.edu/ ~cole/ic 1613. fig8.eps 



Fig. 8. — Synthetic CMDs created assuming a Salpeter IMF and constant star formation rate from 250-1000 Myr. The Hess 
diagram for IC 1613 is overlain. Left panel: model Z = 0.001. Right panel: model Z = 0.004. The red giant branch, older than 
1000 Myr, has not been modelled here so that the locus of young giants and supergiants can be seen more clearly. 



1; secondly, the base of the RSG plume with I <: 22 
and V— I ~ 1. Dohm-Palmer et al. (1997b) presented 
the first unambiguous detection of the blue loop stars 
in a metal-poor system and showed how these stars 
can be used to infer the star formation history of a 
galaxy over the past 1 Gyr. The clearest inference 
from Fig. |8| flows from the termination of the blue 
loop sequence at I « 22; this morphology indicates a 
higher past star-formation rate that declined at some 
point; the age of decline is metallicity dependent: 300 
Myr for Z = 0.004, 400 Myr for Z = 0.001. The ex- 
tent of the blue loop sequence in the models is clearly 
metallicity dependent; we see from Fig. ||] that the Z 
= 0.001 models are bluer than the data, while the Z 
= 0.004 models are too red. The same behavior can 
be seen with respect to the RSG sequence. An inter- 
mediate metallicity, Z rs 0.002, is therefore predicted 
for the 0.1 - 1 Gyr old stellar populations of IC 1613. 
This value is in excellent agreement with the nebular 
abundances of Skillman et al. (1989) and Davidson 
& Kinman (1982). 

We have not directly compared the IC 1613 Hess 
diagram to the WIYN CMD; these data do not reach 
beyond 100 Myr on the main-sequence, and are thus 
unmodelled by the Girardi et al. (1999) isochrone set. 
As in the WFPC2 CMD, the Z = 0.004 models are too 
red when compared to the RSG plume. Even the Z 
= 0.001 model blue loops do not reach far to the blue 
as the blue loops in the WIYN CMD above I m 21.5. 
If the models, and the photometry, can be trusted on 
this point, this would be indicative of a metallicity 
well below Z = 0.001 for the youngest generation of 
stars in IC 1613. However, comparison of the WIYN 



and WFPC2 CMDs shows that crowding errors affect 
even the brighter stars in ground-based images (see, 
e.g., Tolstoy et al. 1998). Additionally, it has been 
shown that the color extent of blue loops in models 
of massive stars is quite sensitive to the adopted in- 
put physics (e.g., Chiosi 1998; Langer 1998); effects 
such as the mass-loss prescription, semiconvection, 
and rotation (e.g., Langer & Maeder 1995) can lead 
to a wide variety of blue loop morphologies at a given 
age and metallicity. 

5.2. Intermediate- Age and Old Stellar Populations 

Baade (1963) and Sandage (1971a) noted the ellip- 
tical sheet of red stars within which the main body 
of IC 1613 is embedded; drawing on Baade's (1944) 
division of the stellar populations of the more lumi- 
nous Local Group galaxies into Populations I and II, 
this sheet has been identified with the Pop II ele- 
ment of IC 1613. Freedman's (1988a) photometry 
resolved this sheet into red giant and aymptotic gi- 
ant branches, to a level of V ~ 23, confirming the 
presence of stellar populations aged at least ml Gyr, 
and perhaps as old as the most ancient Galactic stars, 
10-15 Gyr old. 

Our WIYN color-magnitude diagram extends 
Freedman's (1988a) results; the strongest feature of 
the WIYN CMD (Fig. |) is the classic Pop II red gi- 
ant branch. We detect red giants to a limiting magni- 
tude of I ps 23 (V « 24) across the whole of the 6'.8 x 
6'.8 field of view. Although the center of IC 1613 has 
not seen significant recent star formation and is the 
site of a hole in the H I distribution (Lake & Skillman 
1989), the stars of our WIYN image (Fig. |j) show 
no tendency to avoid this area. This is in agreement 
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with the results of Hodge et al. (1991), who noted 
the differing behaviors of the bright blue stars and 
the faint red stars in this regard. 

Because of the near-degeneracy of the red giant 
branch V— I color to age beyond 2 Gyr, it has not 
previously been possible to infer the age mixture of 
the Baade's sheet population in IC 1613 from CMD 
morphology (Freedman 1988a); to resolve this issue 
requires photometry further down the main-sequence. 
Saha et al. (1992) identified 15 RR Lyrae-type vari- 
able stars in IC 1613; these betray the presence of 
a genuinely ancient stellar population in this galaxy. 
The observed variation in RR Lyrae production ef- 
ficiency among Galactic globular clusters is evidence 
of an unknown ( "second" ) parameter which impacts 
the horizontal branch morphology of ancient stellar 
populations (van den Bergh 1967; Lee et al. 1994). 
This effect limits the power RR Lyrae number counts 
to constrain the precise contribution of ancient stars 
to the IC 1613 RGB; a lower limit of 15% was placed 
on the fraction of red giants older than psIO Gyr. 

5.2.1. The Horizontal Branch 

The presence of a red horizontal branch (RHB) 
population in the WFPC2 VI color-magnitude dia- 
gram can be discerned, in the color range V— I ps 0.6- 
0.8, at I ~ 24.2-24.3. This is consistent with the mean 
g magnitude (Thuan & Gunn 1976) of 24.9 found by 
Saha et al. (1992) for the IC 1613 RR Lyraes. The 
Pop II stars of the central 500 pc of IC 1613 thus seem 
similar to the those in the relatively distant field stud- 
ied by Saha et al. (1992). The theoretical RHBs of 
recent isochrone sets (e.g., Girardi et al. 1999) pre- 
dict a magnitude of I ~ 24.35, for a metal abundance 
of 1/19 solar (Z = 0.001, corresponding to [Fe/H] = 
-1.3) at the distance of IC 1613. The -0.15 mag off- 
set between the observed RHB and model isochrones 
can be accounted for by a difference in [Fe/H] of 
ps -0.5 to -0.7 dex (Caloi, D'Antona & Mazzitelli 
1997), suggesting a metallicity for the RHB popula- 
tion of [Fe/H] ps -1.8 to -2.0. 

The red spur which extends from the main- 
sequence to V-I ps 0.3 at I ps 24.7-24.8 cannot 
be definitively assigned to a specific evolutionary 
state without further analysis. At least two different 
classes of stellar population can create such a feature 
in the color-magnitude diagram. The red spur could 
be post-main-sequence-turnoff stars aged 900 to 1100 
million years with Z = 0.001; alternatively, its color 
and magnitude are consistent with the blue horizontal 
branch (BHB) of a 10-15 billion year old population 
of very low metallicity (Z ps 0.0004). 

The BHB interpretation shows consistency with 
previously known properties of IC 1613: its old aver- 
age age, the inferred metallicity of the red horizontal 



branch stars, and the detection of RR Lyrae variables 
by Saha et al. (1992). It is nigh impossible to pre- 
cisely predict the expected number of BHB stars from 
IC 1613's RGB luminosity function, because the RGB 
is likely contaminated with intermediate-age stars, 
and because the second parameter effect decouples 
horizontal branch morphology from the known phys- 
ical properties of ancient stellar populations. Nev- 
ertheless, we performed a simple analysis to estimate 
the relative star-formation rates required to fully pop- 
ulate the red spur with BHB stars. We used a Monte 
Carlo technique to draw stellar populations from the 
Girardi et al. (1999) isochrones according to the ini- 
tial mass function of Salpeter (1955) from 0.7 - 5 Mq, 
to reproduce the relative numbers of stars seen in the 
red spur and on the main-sequence at the same I mag- 
nitude. From this we infer that the star-formation 
rate 15 Gyr ago was some 2-200 times larger than the 
recent (t <^ 0.9 Gyr) star-formation rate. The limiting 
uncertainty in this calculation is the second parame- 
ter effect (see Saha et al. 1992). The main-sequence 
turnoff associated with this dominant, ancient popu- 
lation would be located at I ps 27.5-28, requiring an 
8-meter class telescope to photometer reliably. 

We used the same analysis technique to evaluate 
the hypothesis that the red spur at I ps 24.7 is the 
subgiant branch produced by a burst of star forma- 
tion ps 1 billion years ago. We divided the stars in 
this region of the CMD (25 > I > 24.3) into two color 
bins: a "main-sequence" bin from —0.4 < V— I < 0.1, 
and a "red spur" bin from 0.1 < V-I < 0.5. The 
number ratio of main-sequence stars to spur stars is 
2.08, which can be reproduced with models in which 
the star- formation rate from 900 to 1100 million years 
ago was a factor of 30 higher than before or since. 

This large increase in star-formation rate would 
have dramatic consequences for the color-magnitude 
diagram of IC 1613, as well as for the galaxy itself. 
A sharp discontinuity in the WFPC2 main-sequence 
luminosity function would be present at I ps 25.1, cor- 
responding to the main-sequence turnoff of the puta- 
tive 1 Gyr burst. Past experience with WFPC2 pho- 
tometry (e.g., Dohm-Palmer et al. 1997a, Gallagher 
et al. 1998) suggests that our photometry is ps50% 
complete at this level. Using all of the WFPC2 pho- 
tometry, we find rising power-law luminosity function 
down to I ps 26 with a sharp decline due to incom- 
pleteness below that level. 

Although increasing photometric errors prevent the 
lower main-sequence luminosity function from ab- 
solutely ruling out a factor of 30 increase in star- 
formation rate a billion years ago, the stars of the 
red clump and giant branch are not compromised by 
these problems. A large 1 billion year-old population 
will produce a very distinctive "vertical red clump" , 
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or pseudo-clump, concentrated in color between 0.7 
<; V-I ^ 0.9 and stretching from M/ « to —1.2 
(Caputo, Castellani & DegFInnocenti 1995; Girardi 
et al. 1998). The model which best reproduces the 
ratio of stars in the red spur to stars on the main- 
sequence predicts a ratio of pseudo-clump stars to 
brighter RGB stars of 6.8, which is inconsistent with 
our observed ratio of 1.5. While IC 1613 shows a 
pseudo-clump, its extension to I & 22 (M/ ss —2.25) 
suggests a period of more or less continuous star for- 
mation from 400-1200 Myr ago. 

A further inconsistency with the 1 Gyr burst sce- 
nario for the red spur is the relative lack of stars 
between V— I i=a 0.3 and 0.6; subgiant stars evolve 
quickly through this region of the CMD, but at a 
smoothly increasing rate. The red edge of the spur 
at V— I 0.3 would imply a sudden acceleration in 
the pace of stellar evolution of low-mass stars which 
is not predicted by theoretical models (e.g., Girardi 
et al. 1999). Contrariwise, the color gap is consistent 
with the location of the RR Lyrae instability strip; 
such gaps are observed in the horizontal branches 
of many systems (e.g., the Carina dwarf spheroidal 
galaxy — Smecker-Hane et al. 1994). 

5.2.2. The Asymptotic Giant Branch 

The strong AGB seen in the WIYN CMD (61 stars 
redder than V-I = 1.8, with I from 19.2-20.0) is 
evidence of a large, intermediate-age stellar popula- 
tion in IC 1613. 17 of these stars are also present in 
the WFPC2 CMD. In the WIYN image, the density 
of AGB stars is highest along the bar of the galaxy, 
and decreases with increasing distance from the main 
body of the galaxy. The AGB density does not in- 
crease towards the areas of current star formation, 
northeast of the bar. 

Cook, Aaronson & Norris (1986) developed a re- 
lation between C/M star ratios and heavy element 
abundances, which was used by Cook & Aaronson 
(1988) to find a metallicity of [Fe/H] < -1.0 for the 
intermediate-age stars in IC 1613. Our WIYN I-band 
AGB luminosity function resembles the C-star lumi- 
nosity function of Cook et al. (1986) above I ~ 20, 
where their sample becomes incomplete. Scaling the 
14 detected carbon stars from Cook et al. to the 
WIYN field of view predicts ~87 carbon stars; be- 
cause the very red stars are somewhat clustered to- 
wards the northwest quadrant of the WIYN image, 
where one of the Cook et al. fields was located, this 
overprediction is not a cause for concern. Conversely, 
it is not unreasonable to expect a large fraction of our 
61 AGB stars to be bona fide carbon stars. Cook et 
al. (1986) found the carbon star fraction among ex- 
tended AGB stars to be 75%-95% in their two fields. 

Mould &: Aaronson (1982) recognized the fact that 



upper AGB stars are a tracer of 2-6 billion year old 
stellar populations, and used the luminosity extent 
of the AGB as an indicator of stellar age. Reid &: 
Mould (1990) examined the AGB population of the 
Small Magellanic Cloud and found that, like IC 1613, 
the SMC has a high fraction of carbon stars among 
the faint (-3 < M bo i < -5.5) AGB stars. From the 
AGB luminosity function, Reid &: Mould concluded 
that the mean age of the SMC is greater than 2-3 
Gyr. 

We perform a similar analysis, adopting the car- 
bon star bolometric correction from Reid & Mould 
(1985), BCi = 1.9 - 0.7(V-I) . The resulting lumi- 
nosity function for our WIYN data is shown in Figure 
|9|. The Mboi-age relation from Table 3 of Mould & 
Aaronson (1982) is marked with open stars and la- 
belled according to age. Incompleteness is increasing 
towards brighter M^oi: the most luminous stars in 
Figure [9| are located at V ~22.5, near the limit of the 
WIYN data. Therefore we may be biased towards the 
detection of somewhat older carbon stars. 

It is clear from Figure | that IC 1613 has formed 
substantial numbers of AGB stars over a long period 
of time. The lack of AGB stars younger than 2 Gyr is 
partially due to the decreasing stellar lifetimes with 
decreasing age. However, a higher star formation rate 
during ancient times than during recent ones is con- 
sistent with the presence of the blue horizontal branch 
stars. Information derived from such an advanced 
phase of stellar evolution is subject to compounded 
uncertainties (see, e.g., Chiosi 1998). Marigo, Girardi 
&: Chiosi (1996), for example, had difficulty reproduc- 
ing the distribution of carbon stars in SMC clusters 
with existing models; even with the latest models, 
Marigo & Girardi (1998) found that the shape of the 
carbon star luminosity function is relatively insensi- 
tive to the form of the star-formation history of a 
galaxy. Therefore we decline to draw quantitative 
conclusions regarding the variation of star-formation 
rate over the past 2-10 Gyr of IC 1613's history. 

5.2.3. The Red Clump 

The relative strengths of the ancient and middle- 
aged stellar populations in dwarf irregular galaxies 
are still very uncertain (e.g., Tolstoy 1998; Da Costa 
1998; Mateo 1998). The strong red clump in Figure ^ 
is an indication of the number of 1-10 Gyr old stars. 
Because the red clump stars share a common core 
mass (M c ~ 0.5 Mq), their helium-burning lifetimes 
do not increase rapidly with increasing age (decreas- 
ing total mass) from ages ~ 1.2-10 Gyr. In contrast, 
the RGB lifetime is a strongly decreasing function of 
stellar mass. The ratio of RC stars to RGB stars in 
a CMD can therefore be used as a rough estimate of 
the mean age of a stellar population (see, e.g., Han 
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Fig. 9. — Bolometric luminosity function of AGB stars in the WIYN field. Results from the carbon star Mboi-age relation 
of Mould & Aaronson (1982) are marked with stars; star formation has been, on average, a continuous process throughout the 
lifetime of IC 1613, although relative star-formation rates cannot be accurately inferred from this figure due to model uncertainties 
and small number statistics. 



et al. 1997, Gallagher et al. 1998, Tolstoy 1998). 

In the WFPC2 CMD of IC 1613, the number ra- 
tio of RC to RGB stars is 2; because photometric 
completeness at the level of the red clump (V, I 
« 24.6, 23.8) is w (80%, 90%), we take our value 
N(RC)/N(RGB) > 2 to be a lower limit. Examina- 
tion of synthetic CMDs therefore provides an upper 
limit to the mean age of the RGB-producing popula- 
tion: t RGB < 7 ±2 Gyr. The ratio of RC to RGB life- 
times has a slight metallicity dependence, such that 
a more metal-rich population is older than a metal- 
poor one at a given N(RC)/N(RGB) (Cole 1999, in 
preparation) . 

We can make a further estimate of the mean age of 
the red clump stars using its mean I magnitude. The 
magnitude of the red clump depends on its metallic- 
ity and, to a lesser extent, its age (e.g., Girardi et al. 
1998; Seidel, Demarque, & Weinberg 1987). To define 
the magnitude of the red clump, we plot the WFPC2, 
I-band luminosity function in Figure 10. The mean 
magnitude of the red clump stars is I = 23.76, and the 
dispersion about the mean is a = 0.33. This mean 
magnitude corresponds to an absolute magnitude Mi 
= -0.56. 

If we adopt the mean RGB metallicity of [Fe/H] 
= —1.4 for the red clump, we can then solve for 
the mean age which produces a clump of the cor- 
rect brightness for that metallicity. This is most con- 
veniently done relative to the Solar neighborhood, 
where Mt(RC) = -0.23 ±0.03 (Stanek & Garnavich 
1998). From Cole (1998), we find that the mean 
age of the IC 1613 red clump is approximately equal 
to that of the Solar neighborhood, t ~ 6 Gyr (e.g., 



Jimenez, Flynn, & Kotoneva 1998). 

Both of these estimates are subject to strong sys- 
tematic error. Contamination of the RC by stars 
in the 800-1200 Myr age range can bias the mea- 
sured Mt(RC) either high or low, depending on the 
exact star formation and chemical enrichment histo- 
ries. Photometric completeness decreases by a small 
amount across the magnitude range of the clump, and 
may induce a slight bias toward younger ages. On a 
more fundamental level, there is still some question 
as to whether the red clump magnitude varies with 
age between 2-10 Gyr as models predict — Udalski 
(1998) argues that the clump magnitude does not 
vary with age in a sample of Magellanic Cloud star 
clusters; however, Girardi (1999) finds that the pre- 
dicted fading with age is not ruled out by Udalski 's 
data. Despite the formiddable observational and the- 
oretical uncertainties in the analysis, we feel that our 
main conclusion — that IC 1613 is dominated neither 
by 10-14 Gyr old stars nor by 1-4 Gyr old stars, but 
contains stars of all ages — is robust. 

We can further use the red clump to constrain the 
metallicities of the stellar populations of IC 1613 at 
various ages. The "pseudo-clump", which overlays 
the red clump but consists of stars just above the 
maximum mass of core-helium flash (1.7-2.0 Mq, 
e.g., Sweigart, Greggio & Renzini 1990) contains stars 
~1-1.2 Gyr old. The color of this feature has a strong 
metallicity-dependence (Girardi et al. 1998); from 
the lower limit of V— I ~0.75 for the red clump we 
find a lower limit of [Fe/H] rs — 1.5 for the 1—1.2 Gyr 
population. 

The fading of the red clump with increasing age and 
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Fig. 10. — RGB I-band luminosity function, highlighting the red clump. The mean magnitude of the clump is I = 23.76 ±0.30. 
The red clump here contains a mixture of post-core-helium flash, "true" clump stars, older than 1.2 Gyr, and more massive stars 
aged 0.8-1.2 Gyr. This combination and growing incompleteness at the faint end complicate the interpretation of the clump. 



metallicity allows us to put an upper limit on [Fe/H] 
for the very old stars in IC 1613, using the magni- 
tude of the horizontal branch and lower red clump. 
For ages older «6 Gyr, the Z = 0.004 ([Fe/H] = -1.3) 
isochrones from Girardi et al. (1999) become too faint 
to match the observed red clump in IC 1613. This 
is essentially the same effect that we have seen for 
the red horizontal branch stars, which have [Fe/H] ~ 
-1.9. 

6. DISCUSSION 

The analysis of stellar populations allows us to con- 
struct a rough age-metallicity relation for IC 1613. 
Previous determinations of IC 1613's global metal- 
licity have mainly been from nebular spectroscopy 
(e.g., Skillman et al. 1989). These measures sample 
only the youngest component of IC 1613, and have es- 
tablished the current metal abundance of IC 1613 at 
[Fe/H] « -1. Freedman (1988a) assessed the RGB 
abundance, sampling the stars from 2-14 Gyr old; 
Cook & Aaronson (1988) judged the upper limit of 
[Fe/H] to be < — 1 for the 2-6 Gyr population, using 
the ratio of C to M stars along upper AGB. 

In this paper, we have estimated the metal abun- 
dance of the blue and red supergiants to be [Fe/H] 
~ —1, providing an abundance indicator for stars in 
the 10-500 Myr range. We reappraised the metallic- 
ity of the RGB, finding [Fe/H] = -1.4 ±0.3. Using 
the morphology of the red clump, we provide an up- 
per limit of [Fe/H] <^ —1.3 for the stars older than 
~ 6 Gyr. And from the magnitude of the horizon- 
tal branch stars, we infer a metallicity for the old- 
est stars in IC 1613 of [Fe/H] ~ -1.9. We have as- 



age-metallicity relation for IC 1613. Overplotted is 
an age-metallicity relation for the Small Magellanic 
Cloud, from Pagel & Tautvaisiene (1998). Despite 
the large uncertainties, it appears that the chemical 
evolution of IC 1613 has taken a similar form to that 
of the SMC, although IC 1613 may be slightly more 
metal-poor. Figure 11 quantifies to an extent the 



sembled these data into Figure 11, our preliminary 



"population box" diagram given in Hodge (1989), al- 
though we find a stronger trend of metallicity with 
age than the population box constructed from older 
data by Grebel (1998). 

Where do IC 1613's stellar populations place it 
within the context of the Local Group dwarf galax- 
ies? Many Local Group dwarfs have now been stud- 
ied with WFPC2 to a comparable depth as IC 1613, 
and so we can begin to use direct intercomparisons 
to test the accuracy of model-dependent interpreta- 
tions of CMDs in terms of galaxy evolution. In Figure 
|T^ , we plot the grayscale Hess diagrams of four Local 
Group dwarf galaxies, and overplot the Hess diagram 
of IC 1613 for comparison. The basic parameters of 
each galaxy, taken from Mateo (1998), are given in 
Table 2. The quantities A(m— M) and AEb-v de- 
note the offsets applied to the data for each galaxy 
in order to compare the Hess diagrams to IC 1613. 
Each diagram is constructed from WFPC2 F555W 
and F814W data, and plotted with a logarithmic 
transfer function. 



6.1. NGC U7 

NGC 147 is a dwarf elliptical satellite of M31; the 
WFPC2 data are from Han et al. (1997). Its RGB 
extends further to the red than IC 1613's, due to the 
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Fig. 11. — Age-metallicity relation for IC 1613. Each region of the figure is labelled with the CMD feature that constrains the 
metallicity to lie within the shaded region. Abbreviations are as given in the text; WR denotes the Wolf-Rayet star; HII denotes 
the H II regions. The dotted line shows the mean age-metallicity relation for the SMC. 



Table 2 
Four Local Group Dwarfs 



Galaxy 


A(m-M) 


AEb-v 


M v 


[Fe/H] 


Mhi/Lb 


IC 1613 






-15.2 


-1.4 


0.81 


NGC 147 a 


0.05 


0.15 


-15.5 


-1.1 


<0.001 


Pegasus b 


0.5 


0.12 


-12.9 


-1.3 


0.44 


Leo A c 


0.0 


0.0 


-11.4 


-1.7 


1.6 



a Han et al. 1997. 
b Gallagher et al. 1998. 
c Tolstoy et al. 1998. 
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Figure 12 

can be obtained from 
the following URL: 

http:/ /www. astro, wisc.edu/~cole/ic 16 13. fig 12.eps 



Fig. 12. — Greyscale Hess diagrams for 4 Local Group dwarf galaxies of similar luminosity: IC 1613, NGC 147, Pegasus, and 
Leo A. The IC 1613 logarithmic Hess diagram is contoured atop each panel to highlight the similarities and differences between 
them. 
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higher mean abundance (and/or age) of NGC 147. 
Concordantly, the NGC 147 red clump lies «0.1 mag 
below IC 1613's. The upper main-sequence is unpop- 
ulated in NGC 147, a sign that star-formation ended 
at least «1 Gyr ago (Han et at 1997). A blue hor- 
izontal branch population is prominent in NGC 147, 
sharing its locus in the Hess diagram with that of 
IC 1613. In NGC 147, star-formation went nearly to 
completion early in its history, enriching subsequent 
generations of stars, which now populate the RGB. 

6.2. Pegasus 

As mentioned above, Pegasus is a transition galaxy 
between the dwarf irregular and dwarf spheroidal 
classes; it is also one of IC 1613's nearest neighbors 
within the Local Group. Though less massive and 
less luminous than IC 1613, Pegasus is of compara- 
ble metal abundance to IC 1613, and has a similar 
RGB and main-sequence morphology. However, Pe- 
gasus has a larger AGB population than IC 1613, as 
would be expected from its enhanced period of star- 
formation a few Gyr ago (Gallagher et at 1998). The 
recent star-formation rate in the WFPC2 field, M ~ 
1600 Mq Myr -1 kpc -2 , is similar to that derived 
by Gallagher et at (1998) for Pegasus, M = 1350- 
2000 M Myr" 1 kpc -2 ; however the WFPC2 image 
of IC 1613 excludes the most active regions of star 
formation. Pegasus and IC 1613 may be examples of 
galaxies which undergo occasional strong episodes of 
star-formation lasting up 0.5-1 Gyr, at intervals of 
one to a few Gyr. The differences in their recent (<3 
Gyr) star-formation histories would then be due to 
a difference in phase in their cycles of quiescent vs. 
elevated star-formation. 

6.3. Leo A 

Leo A provides a stark contrast to the other galax- 
ies in this sample. The color-magnitude data are 
from Tolstoy et at (1998). The smallest, most iso- 
lated galaxy in the sample, the CMD of Leo A shows 
the unusual feature of a dominant blue supergiant 
branch, evidence for strong recent star formation. 
Leo A's giant branch lies to the blue of IC 1613's, in- 
dicating its lower metallicity (and/or age). The blue, 
vertically extended red clump in Leo A suggests that 
this galaxy is dominated by an ~ 1 Gyr population. 
Deep WIYN & MMT images of Leo A have failed to 
reveal an extended Baade's sheet population of red gi- 
ants that would be indicative of significant early star 
formation (A. Saha, private communication, Tolstoy 
et at 1998). The low metallicity and high gas frac- 
tion of Leo A are all consistent with a galaxy that has 
not had a major star formation episode until quite 
recently. Leo A is even more isolated than IC 1613 
(its nearest neighbors are the Leo I and Leo II dwarf 



spheroidals, more than 450 kpc away), and it must 
stand as a counterexample to the Tucana dwarf, a 
purely old system (Da Costa 1998), in evolutionary 
studies of small, isolated galaxies. 

7. SUMMARY 

We have used WFPC2 to obtain deep images of the 
dwarf irregular galaxy IC 1613. We analyze the re- 
sulting (I, V— I) color-magnitude diagrams to study 
the stellar populations at the center of this object. 
The strongest feature of the CMD is the red giant 
branch and superposed red clump. We find a mean 
metallicity of [Fe/H] = —1.4 ±0.3 from the color of 
the RGB, consistent but slightly lower than previ- 
ous results (Freedman 1988). Upper main-sequence 
stars are also present, indicating a continuous period 
of star-formation during the past ~350 Myr at a level 
of ~ 3.5 x 10 -4 M yr" 1 per WFPC2 field, corre- 
sponding to 1600 Mq Myr -1 kpc -2 . From a compar- 
ison of the supergiant stars to theoretical stellar evo- 
lution models, the mean chemical abundance during 
the past 800 Myr has been approximately one-tenth 
solar, consistent with previous measurements of H II 
region abundances. 

In our WIYN image, we find 61 candidate carbon 
stars which populate an extended upper asymptotic 
giant branch and trace intermediate-age (2-10 Gyr) 
stellar populations. The AGB luminosity function 
indicates broadly continuous star formation from at 
least 4-10 Gyr ago; a possible decline since then can 
be neither confirmed nor ruled out with the present 
data. 

We have detected for the first time a blue horizontal 
branch population in IC 1613. This very old popula- 
tion is not indicative of the IC 1613 RGB in general: 
from an analysis of the red clump stars, we find a 
mean age for the RGB of t <J 7 ±2 Gyr. Systematic 
effects are unlikely to push this value upwards by the 
factor of ~2 necessary to enforce coevality with the 
blue horizontal branch. 

We present a preliminary age-metallicity relation 
for IC 1613. While crude, this relation shows the 
same form as that of the Small Magellanic Cloud. We 
will use the age-metallicity information in a detailed 
modelling of the IC 1613 CMD in a subsequent pa- 
per. A comparison to WFPC2 CMD's of the galaxies 
Leo A, NGC 147, and Pegasus shows that the stellar 
populations of IC 1613 are very similar to those of 
Pegasus. Pegasus has an excess of bright AGB stars 
indicating elevated star formation a few Gyr in the 
past, while IC 1613 shows a higher level of current 
star formation; but both galaxies seem to have expe- 
rienced a relatively constant, low level of star forma- 
tion throughout their histories. Relatively isolated 
Local Group dwarf galaxies display a considerable 
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range in their star-formation histories, from objects 
where star formation was completed several Gyr in 
the past (e.g., Tucana — Lavery et al. 1996; possibly 
NGC 147- Han et al. 1997), to young-star dominated 
systems such as Leo A (Tolstoy et al. 1998). 
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